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ABSTRACT 

Approximate solutions of the equations of motion governing laminar 
incompressible fluid flow through a cylindrical channel with a porous 
wall are derived. The invalidity of an approximation in the solution 
of these equations under certain circumstances is pointed out, and the 
results of a numerical integration in the region where the approximation 
is invalid are indicated. A description is given of an experiment 
to verify the calculations, and some interesting results are noted.. 
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INTRODUCTION 


The flow of fluids through channels with porous walls has been of 
interest in the last decade or so because of applications to transpira¬ 
tion cooling, gaseous diffusion, and boundary-layer control. Several 
investigators V ~ ' have contributed much to the theoretical study of 
this type of flow, but little has been done experimentally. 

In general, the theoretical studies have been limited to laminar, 
incompressible flow. Most of the work has been done on the problem of 
a rectangular channel, and some interesting techniques have been developed 
for the solution of this problem. However, in this case it is rather 
difficult to check the solutions by a reasonably precise experiment. 

We chose the more practical experimental problem of flow through 
a finite porous cylindrical tube, accepting the theoretical difficulties 
which are presented. 
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THEORY 


Laminar, incompressible fluid flow through a channel with porous 
walls is governed by the specialized Navier-Stokes equations, the equa¬ 
tion of continuity, and appropriate boundary conditions. 


p(v*V)v - T](V.V)V + Vp = 0, 


V*v = 0. 

p is the density of the fluid and r) is its viscosity. 

For a cylindrical channel with a uniformly porous wall the problem 
reduces to a two-dimensional one described by the following equations: 

d 2 t 


, du du\ /0 2 u , 1 du 

p ( u §? + ' , ' ( ri + 


, 1 Su _ ^ t ^U \ _|_ _ Q 


dF 


dz* 


, dv dv-, /d 2 v , 1 dv , djvq 

P< U 5F + - try + r 3? + 


dp 


+ h M ’ °- 

u = u(r,.z) is the radial component of the velocity, and v = v(r,z) is 
the axial, or z-component of the velocity. 

The boundary conditions arise from the following physical considera¬ 
tions. Because of symmetry there should be no radial velocity on the 
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axis of the tube, i.e. u(o,z)=0. The axial velocity profile must be 

=0. At the wall of the 

r = o 

channel ( r = a) the velocity should have only a radial component. Thus, 
v(a,z) =0. If the flow density is the same through all portions of the 

porous wall, u(a,z) = u . 

The analysis is simplified by the transformations R = r/a, Z = z / a > 
P = p/pu Q 2 , V = - v/u Q , U = - u/u Q , and e = - T)/pu 0 a, which give the 
following dimensionless equations: 


Sv 

flat at the center of the tube, so ^ 


<°§ +v l> - e( iS 


+ 1 dU _ + o_u) + 

R oR 



^(EU) + ^(RV) . 0 . 

The boundary conditions become V(l,Z) = 0, U(l,z) = -1, U(0,Z) = 0, and 

^ = 0. 

R = 0 

A consequence of the assumption that the flow density at the porous 
wall is uniform is that the Z-component of the velocity may be written 
as Z times some function of R, with Z = 0 defined by the point at which 
V(R,Z) = 0. By using this result and introducing a stream function t 

dt pn _ ^ 

that identically satisfies the continuity equation, RV = ^ R , - KU - gg, 
it is found that the expression for the pressure has the form 
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p = p (R) + |yZ 2 , where 7 is a constant for a particular value of the 
parameter e. P (r) may he determined explicitly once the differential 
equations have been integrated. 

A particular stream function which satisfies the continuity equation 
identically and. also fits the boundary conditions is f = BZf(R 2 ), where 
B is a constant. The argument of the function f is chosen to he R 2 
rather than R for convenience in analysis. 

The substitution X = R 2 gives ty = BZf(x), U = - 1 ^ (x), and 

VX 

V = 2BZf'(X), where the prime denotes differentiation with respect to 
the argument of the function, which in this case is X. The velocity 
profile is given by V(X). 

With the above expressions for U and V the differential equations 
become 



2B 2 ff1 

nr 

_ b2 f 2 + beB nr f” + 2 /T" P ' = 0 , 

x 3/ 2 

(1) 

and 

B 2 (ff ' 

' - f ,2 ) + 2eB(f " + Xf' ") = J • 

( 2 ) 


If f is known explicitly, then P q is determined within a constant of 
integration. 

P q = - B 2 dX + f J ^ clX - 2eBf' + K. 

Equation ( 2 ) may be re-formulated for attack in two limiting cases. 
In the limit of large e 

|(f> ' + Xf " ') - -(f ,2 - ff ") = - 2 — = b. (3) 

B 6 4eB 2 
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This may be regarded as the equation governing the flow of fluids with 
high viscosities, or as the equation governing very slow flows. For 
high flows or small viscosities the limit of small e gives 

f < 2 - ff " - §£(f " + Xf *' ’) = - — = c. (4) 

B kB S 


Equations ( 3 ) and (4) are subject to the boundary conditions f'(l) 0, 

V X 


Equation ( 3 ) may be expanded around (l/e) = 0 to obtain a pover 
series solution in (l/e). Since 7 may be a function of e^ d must also 
be expanded around (l/e) = 0. For greater clarity the solution of 
equation ( 3 ) will he denoted by g(X), while the solution of equation (if) 
will remain f(x). 


8 = 6 0 + hx + ^ 2 g 2 + ‘ 

d = d + id + —d + . 

O €1 2 2 

e 


> 


For this series expansion the boundary conditions become 
g 0 (D - | , g n (l) 

g n ,(l) = °' n - °' 11m = 0 n ^ 0 

x "° rr~ ’ 

iS /x” g "(x) = 0 , n * 0 . 

Equation ( 3 ) has been solved to order l/e 3 > 'but the solution is 
too unwieldy to present here. To order l/e 2 
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Bg(X) = - (X 2 



1/X? xf. 

+ £ l 36 “ 6 



1 ,x 6 X 5 X 4 

2 5400 “ 360 72 


11 X 3 . 19X 2 83 X 

216 + 270 " 2700 


88 1 187 1 

and 7 = - l6e - 12 + ^ € - 315 


) 


In a like manner equation (4) may be expanded around e 0. 


f = f + ef + e 2 f 
0 1 2 

c = c + ec + e 2 c + . . . . 
o 1 2 

A similar set of boundary conditions holds. 

To zero order in e equation (1) has the solutions 


.... Recall that 

X = r 2 = r 2 /a 2 . These functions give a corresponding 7 - " ( 2n + l) tt 2 . 
The fact that there are an infinity of solutions may be tied in with the 
boundary condition at the exit end of the tube. The only solution that 

has been seen experimentally is that for which n = 0. 

The solution of equation 00 has been carried out to order one in 
e ^ 2 ) p or the flows investigated the contribution of the first order 
term is small. At the lowest flow (largest e) which gives useful data 
this term gives less than a 1* correction, and its effect.decreases with 

decreasing e. 

It should be pointed out that the zero order small e approximation 
amounts to neglecting the highest order derivative, and this term may or 


f = W 

o B 


(2n + 1 )ttX| 


n = 0, 1, 
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may not be important. Numerical integration of an analogous problem 

shows that this procedure is valid when the wall is a source (u^ < 0 )j 

but t hat the highest order derivative dominates near the wall when there 

is suction (u > 0). In the limit of very high flows for wall suction 
o 

the solution is f = X/B. This solution does satisfy the boundary 
conditions and gives a value of 7 of -4. 
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APPARATUS 


Two separate sets of experiments were performed in an effort to check 
the theoretical predictions. The first set yielded qualitative agreement 
with the theory and information necessary for the design of the second, 
more quantitative set of experiments. The porous system used in the 
early experiments, shown in Fig. 1, consisted of an alundum tube which 
was jacketed by a brass cylinder. Six pressure probes evenly spaced 
along the length of the alundum tube at its inside wall were used to 
measure the pressure drop down the channel. Differential measurements 
of pressure against the pressure at the top of the tube on its axis were 
read on a water manometer system. 

This system qualitatively verified the parabolic pressure profile 
(p = -gyZ 2 ), but quantitative verification (7 = - tt 2 ) was not possible. 
Also, attempts to measure the velocity profile proved futile because of 
the sma .11 size of the channel. A larger porous tube that more nearly 
satisfied the theoretical assumptions had to be found. 

The main criteria in selecting a porous tube were that the material 
must have uniform porosity and that it present a suitably large impedance 
to the flow. Graphite, electroplated metal screen, pressed metal screen. 



PRESSURE 

GAUGES 



EXIT 


Fig. 1 - Schematic of Experimental Setup 




sintered metal powder, filter paper, and several forms of alundum 
refractory were tested for these characteristics. Alundum was chosen 
again for the final experiments "because, of the readily available mate¬ 
rials, it most nearly satisfied the requirements. A thicker-walled tube 
was selected' so the impedance would be large enough. 

The tube in the modified apparatus was 18 inches long and had a 
1-1 inch bore and a £ inch wall. Pressure probes to measure y were 
installed only at the top and bottom of the channel. Measurements of 
the velocity profile were made by a pair of differential impact tubes 
which faced into the gas stream. These are shown in Fig. 2. One of the 
tubes was fixed on the channel axis and the other was adjustable radially. 
The impact tubes were mounted in a holder to form a unit which was 
movable along the length of the channel. Accurate knowledge of the 
radial position of the movable probe was lost if the unit was moved 

away from the bottom of the tube. 

In both experiments the gas flow was up from the bottom of the 
vertical unit, through the porous wall, and back to the bottom where it 
was discharged to the atmosphere. 
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Fig. 2 - Velocity Profile Impact Tubes 










RESULTS 


The limiting value of the parameter 7 and the magnitude and shape 
of the velocity profile are the only things that have been successfully 
checked experimentally, and these were found only in the case of small 
e (large flows). Experimental error was too great to allow a determina¬ 
tion of 7 as a function of e, and the sensitivity of the measurements 
could not be made high enough to give meaningful results at low flows. 

The theory shows that when the wall is a source (u Q < 0) 7 should 
have a limiting value -of - tt 2 , and the velocity profile should be given 

*y 

V = 2BZf' (X) = ttZ cos jfxj . 

The axial velocity as a function of the position variables is then 


v = 


ITU z 

__ 0 _ 

a 


cos 


TT 




7 and v must now be expressed in terms of measurable quantities. 

Note that 7 = 2AP/Z 2 , where AP = Ap/pu 2 and Z = z/a. Therefore 


7 = 


2 a g Ap 

pu 2 z 2 
o 


Q, in gm/sec. 


u may be expressed in ‘berms of the total measured flow, 
o 

through a conservation of mass relationship. If L is the 
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length of the porous tube, the mass of fluid into the vail each second 
betveen z = 0 and arbitrary z is Qz/L. The quantity out of the wall 
every second is JA, where J is the flow density, -pu Q , and A is the 
cylindrical area of the wall section considered. Thus 
Qz/L = JA = -pu Q 2iraz, 


and 



Q 2 

4iT 2 a 2 L 2 p 


This gives the expression 7 = 8 rr 2 a 4 L 2 pAp/Q 2 z 2 . In the experiment Ap is 
the pressure diffsrence between the points z = L and z = 0. Therefore 


7 = 87T 2 a 4 p— . 

Q 2 


(5) 


Ap, and hence 7 , can be shown to be a negative number by considering 
Bernoulli's Law along a streamline. 

The differential impact tubes measure the difference 

^r = P o " P r = ^ (v o 2 “ V } > 

where v and v are the fluid velocities in the axial direction on the 
o r 

axis and at an arbitrary radial position, respectively. Since it is incon 

venient to measure v for each flow, and since the velocity at the wall, 

o 

v , is small, v is eliminated in favor of v . 
w o v 


AP = P 

w o 


P = |p(v 2 - V 


w 


w 


! ) 


Thus v = / — /"AP - AP - P . In the experiment P was found to be 
*rvp v wrw w 

negligibly small, so the velocity profile is expressible in terms of two 
differential pressure measurements: 
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(6) 



The subscript on the AP's indicates the position of the movable impact 
tube. 

Comparison of theory and experiment for the velocity profile is 
facilitated if the shape and amplitude of the profile are given sep¬ 
arately. For the shape, equation (6) is normalized so the data may be 
directly compared with cos | . The theoretical amplitude is 

^ M 

V = - - = - } ^ ' ' 

°t a 2pa 2 

and the experimental amplitude is given by the normalizing factor 

v - [! TET~ . 

o e V P w 

When the wall is a sink (u q > 0) the boundary-layer type of flow 
that results has a flat velocity profile (except near the wall). 7 W 
be calculated knowing that Ap/pv 2 = - i for this type of profile. Using 
a mass balance equation for the axial direction, Qz/L = JA = pv/ra 2 . 

= l = - \ • 

pv 2 Q 2 Z 2 

It is seen that 7 = -4 for u q > 0. 

All quantities in equations ( 5 ), ( 6 ), and ( 7 ) are in cgs units. 

However, since it is more convenient to measure Ap and the AP's in units 

of mm of water, a correction factor of 98 ^^f^ter mUSt be introduced 
This, along with the 19^ mm radius of the porous tube, gives the results 


20 


(5a) 


7 = 1.063 x 10 5 p^ , 

Q 2 

- ^ . 




V 


r = £ii2WQf . 

0 p L 


( 6 a) 


(7a) 


When the wall is a source^ p is the density of the gas used (nitrogen) 
at room temperature and Los Alamos atmospheric pressure. This density 
is 0.970 x 10" gm/cm 3 . Thus 

(5b) 


Ap 

7 = 103 — , 


\ = 449.5yAP w - AP , 


v = 139 Q T - . 
°t L 


( 6 b) 

(7b) 


The Ap and the APIs are in mm of water. 

Only 7 may be checked quantitatively when the wall is a sink, and 
in this case an additional measurement to find the density of the gas 
must be made. In terms of the absolute operating pressure inside the 
porous tube, II, psi, 

7 = 9.0511 — . (5°) 

Q 2 

The experimental values of 7 were close to the theoretically 
predicted values for both flow directions. It was found that when the 
wall is a sink, 7 = - 3.88 + 0.08 (- 4.00 was predicted from the theory), 

and when the wall is a source, 7 = -9.88 + 0.11 (- 9.87 was predicted). 
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These results are the averages of 15 and 14 determinations, respectively, 
with standard deviations, over a range of flows from 5.84 to 47.7 gm/sec 
of nitrogen. The experiments were done at a pressure of 11.6 psia 

and a temperature of about 25 °C. 

The amplitude of the velocity profile for u q > 0 could not be 
found experimentally, but the flat character of this profile was verified 
within the sensitivity of our pressure measuring system (a water manom¬ 
eter). 

The velocity profile data for u^ ■s 0 are presented in Table 1, with 
a "typical" set superimposed on the theoretical profile in Fig. 3. 
Experimental precision on the steep slope portion of the curve was not 
as high as on the flatter portions because small errors in knowledge of 
the position of the movable probe change the value of the cosine con¬ 
siderably. This introduced considerable experimental error when the 
probes were moved away from the bottom of the tube (z/L = l). 

In Table 1 the velocities given are in meters/sec. The discrepancies 
between the theoretical and experimental velocities may be explained at 
least in part by the finite size of the impact tubes used to measure the 
velocity profiles. Thus it was not AP^. that could be measured, but some 
AP, at a point not on the wall, which was too low. 

Considerable effort was expended in trying to force one of the 
"higher-order" modes of flow to appear, for which n > 0, but these 
experiments failed. The velocity profile for the n = 1 mode is shown 
in Fig. 4. The general method of attack on this problem was to try to 
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Fig. 3 - Velocity Profile for the n = 0 Mode 




VELOCITY PROFILE DATA ANALYSIS 
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Fig. 4 - Theoretical Velocity Profile for the n = 1 Mode 




control the exit boundary conditions so they would be favorable for the 
establishment of the desired flow scheme. The shape of the stream 
function was calculated and a plug which fit this function was machined 
and inserted into the gas stream. Since the experiments failed it has 
not been determined whether or not these "higher-order" solutions have 
physical significance. 

It is interesting to note that the equations of motion that are 
solved are those that govern laminar, incompressible fluid flow, while 
the experiments were performed at exit Reynolds numbers high enough (up 
to about 100,000) that turbulence should have been well established, 
and with nitrogen, a compressible gas, as a working fluid. This indicates 
that the transition from laminar to turbulent flow either takes place 
only at even higher Reynolds numbers, or has little effect on the velocity 
profile. 

Compressibility effects were quite small under the conditions used 
in our experiments. The maximum Mach numbers attained were 0.16. 
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